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ABSTRACT 

We present Herschel PACS 100 and 160 iim observations of the solar-type stars a Men, HD 88230 and HD 210277, which form part of the FGK 
stars sample of the Herschel Open Time Key Programme (OTKP) DUNES (DU si around A^£^arby 5" tars). Our observations show small infrared 
excesses at 160 fim for all three stars. HD 210277 also shows a small excess at 100 jum, while the 100 //m fluxes of a Men and HD 88230 agree 
with the stellar photospheric predictions. We attribute these infrared excesses to a new class of cold, faint debris discs, a Men and HD 88230 
are spatially resolved in the PACS 160 //m images, while HD 210277 is point-like at that wavelength. The projected linear sizes of the extended 
emission lie in the range from ~ 115 to < 250 AU. The estimated black body temperatures from the 100 and 160 jum fluxes are < 22 K, while 
the fractional luminosity of the cold dust is L^ustlL^ ~ 10~^, close to the luminosity of the Solar-System's Kuiper belt. These debris discs are the 
coldest and faintest discs discovered so far around mature stars and cannot easily be explained by invoking "classical" debris disc models. 

Key words. - Stars: planetary systems: planetary discs -Stars: planetary systems: formation - Stars: individual: a Men (HIP 29271) - Stars: 
individual: HD 88230 (HIP 49908) - Stars: individual: HD 210277 (HIP 109378) 



1. Introduction 

Debris discs are tenuous structures associated with main se- 
quence stars formed by second generation dust, which has re- 
sulted from the collisions of solid bodies continuously supplying 
the circumstellar environment with small dust particles. This for- 
mation sequence is inferred from the lifetime of the dust grains 
against destructive collisions, Poynting-Robertson drag and ra- 
diation pressure, which is much shorter than the ages of the host 
stars. These discs are visible in reflected light at optical wave- 
lengths and in thermal radiation at mid-/far-IR and submillimeter 
wavelengths (Aumann et al., 1984; Backman & Paresce, 1993). 
General disc characteristics are grain black body temperatures of 
~ 50 - 100 K, fractional luminosities / > Ldust/Li, ~ 10"^, and 
radii from less than 10 AU to several times 100 (e.g. Absil et al. 
2006, Su et al. 2005, Trilling et al. 2008). Debris discs are con- 
sidered analogues of the Solar System asteroid and Kuiper belts, 
although their luminosities are usually more than 100 times the 
Kuiper belt level of Ldust/Lsun ~ 10"^ - 10"^ (Stern 1996, 
Vitense et al. 2010). The sensitivity of the 3.5 m Herschel far- 
infrared space telescope (Pilbratt et al. 2010) with its instrument 
PACS (Poglitsch et al. 2010) oflfers the possibility of characteris- 
ing colder (~ 30 K) and fainter (L^^^^/L^ few times 10"^) debris 
discs with spatial resolution ~ 60 AU (FWHM) at 10 pc, i.e., 
true extra- solar Kuiper belts. 

DUNES is a Herschel OTKP designed to detect and char- 
acterise extra- solar analogues to the Kuiper belt around main 
sequence FGK nearby stars (Eiroa et al. 2010). In this letter 
we present the results for three stars from the DUNES sam- 
ple: a Men, HD 88230, and HD 210277 as clear examples of 

Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Table 1. Stellar properties. 



Star 


a Men 


HD 88230 


HD 210277 


Sp. Type 


G5V-G7V 


K6V-M0V 


GOV, G7V-G9V 


Distance (pc) 


10.2 


4.9 


21.6 


L*(Lo) 


0.85 


0.15 


1.10 


T.//(K) 


5590 


3850 


5540 


Age (Gyrs) 


5.5 


6.6 


6.9 



the advantages oflTered by Herschel observations: they trace a 
new class of cold, T < 22 K, spatially resolved debris discs with 
very low fractional luminosities. These discs have remained un- 
observed by previous far-IR and submillimeter studies. Table 1 
gives some properties of the stars. Ages are based on the log 
R^^ activity index and have an uncertainty of 60% (Mamajek 
& Hillenbrand 2008). HD 210277 hosts a Jupiter-like planet 
(Marcy et al. 1999). This star and a Men have faint stellar com- 
panions, but neither the measurements nor the photospheric pre- 
dictions are aflTected by them. Eiroa et al. in prep will present a 
full discussion of the stars and the general results of the DUNES 
survey. 

2. Observations and data reduction 

a Men, HD 88230 and HD 210277 were observed with PACS 
100/160 pm in scan map mode. For each star, two scans at po- 
sition angles 70° and 110° were carried out, each scan consist- 
ing of 10 legs with separation of 4'', length of 3' and medium 
speed of 20'7s. Table 2 gives the scan identification numbers 
(Obs. ID) and the total duration of the observations (OT). a 
Men was observed twice in order to increase the S/N ratio. Data 
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Table 2. Log of the PACS 100 yum and 160 yum observations. 



Star 


HIP 


Obs. ID 


OT (sec) 


a Men 


29271 


1342203297/8 


1116 


a Men 




1342216043/4 


2244 


HD 88230 


49908 


1342210610/1 


4500 


HD 210277 


109378 


1342211126/7 


4500 



reduction was made using the Herschel Interactive Processing 
Environment (HIPE) version 7.2. The individual scans were mo- 
saiced to produce the final image at each band. To check the con- 
sistency of the reduction and analysis (particularly regarding the 
eff'ect of correlated noise), mosaics were produced at both the na- 
tive 3V2 for 100 [im (green) and 6! A for 160 yum (red), as well as 
super-sampled pixel scales 1.0'' (green) and 2.0'' (red), the latter 
being the default pixel size in HIPE for this type of data. The 
images with the native pixel size scales avoid, at least partly, the 
correlated noise in the PACS image^(see also Fruchter & Hook, 
2002). A high-pass filter was used to remove large scale back- 
ground emission from the images, with filter widths of 15" and 
25" in the green and red channels. To prevent the removal of any 
extended structure near the stars, regions where the sky bright- 
ness exceeded a threshold value in the image defined by the stan- 
dard deviation of all the positive pixels {S > 10~^ Jy / pixel) 
were masked from that process. Absolute flux calibration un- 
certainties are -3% and -5% for the green and red bands (see 
technical note below). 

3. Results 

Table [3] gives the J2000.0 optical equatorial coordinates of the 
stars as well as their 100 jim peak positions corrected for the 
proper motions of the stars (van Leeuwen, 2007). Ofl'sets be- 
tween the optical and the PACS positions (column 4 of Table [3| 
are within ~ l.Scr the Herschel pointing accuracy of 2'.'4 in this 
observing mode (Sanchez-Portal, private communication). 

Fig. 1 shows the 100 and 160 yum images and isocontour plots 
of the stars, a Men and HD 88230 are point-like at 100 yum, 
FWHM ~ 6'.' 3 X 6 '.'5, while both stars are resolved at 160 yum, 
with angular sizes/position angles ~ 18"/52° and ~ 23'.' 5/45°, re- 
spectively. The extended emission appears clearly asymmetric in 
HD 88230, with the star located at the North-Eastern side. None 
of the objects are resolved in the direction orthogonal to the ex- 
tended emission. HD 210277 is unresolved at both wavelengths. 

PACS fluxes (Table [3]) have been estimated using circu- 
lar and rectangular aperture photometry, taking special care to 
choose the reference background region due to the presence of 
field objects. Specifically, measurements of HD 210277 take into 
account the presence of the bright object located at ~16" North- 
East from the star; in the case of this star we have, in addition, 
carried out PSF photometry using the DAOPHOT software pack- 
age. The PSF photometry fluxes are F(IOO) = 8.4+0.3 mJy and 
F(160) = 14.3+0.4 mJy (errors are those of the PSF fits), which 
are consistent with the aperture photometry estimates. Errors 
have been estimated using a variety of methods, including cir- 
cular and rectangular boxes at diff'erent positions in the nearby 
fields. Sky sizes for the errors estimates are equal to the area of 
the aperture used for the photometry for all three stars and both 
bands; in particular the sizes of the extended emission around a 
Men and HD 88230. 

PACS fluxes have been compared to predicted stellar photo- 
spheric fluxes (Table [3]) using Gaia/PHOENIX models (Brott & 
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Hauschildt, 2005), with the stellar parameters as given in Eiroa 
et al. (in prep). Fig. 1 shows the spectral energy distributions 
(SEDs) of the stars, where PACS fluxes are plotted together with 
optical, near-IR, IRAS, AKARI, and Spitzer/MIFS and IRS data; 
in addition, the best;^-square photospheric fit is shown. To assess 
the presence of an excess at 100 and/or at 160 yum we require that 
the observed fluxes, i^pAcs, exceed by at least 3cr the predicted 
photospheric fluxes, (xa = (^pacs -^*)/<^/i)- No excesses are 
detected at 100 yum for a Men and HD 88230, while it is seen in 
HD 210277. All three stars do show 160 yum excesses. We also 
note that the SED slopes from 100 to 160 jim are a = - 0.5+0.7 (a 
Men), a = - 0.7+0.6 (HD 88230), and a = 0.8+0.8 (HD 210277), 
which clearly difl'er from the expected Rayleigh- Jeans behaviour 
(a = - 2.0) of a stellar photosphere in this wavelength regime. 

3.1. Comments on the PACS images 

All three fields show ~ 2-3 red sources per square arcmin, which 
are likely to be background galaxies. In particular, there is a faint 
peak at -13" towards the SW from HD 88230 and a bright one at 
-16" towards the NE from HD 210277. We have consulted the 
NASA/IPAC extragalactic database to search for counterparts 
without finding any association. Although we cannot firmly ex- 
clude a coincidental alignment or contamination of a background 
source(s) in the line of sight of our stars, we think that it is un- 
likely for the Herschel sources presented here due to the close 
correlation between the optical and Herschel positions, and the 
photospheric predictions and the estimated Herschel 100 jim 
fluxes. In fact, following in a first approach the source counts by 
Berta et al. (2010), the average density of extragalactic sources 
with fluxes ~ 6-7 mJy and ~ 12 mJy - i.e., the measured excesses 
at 160 fim from a Men/HD 88230 and HD 210277, respectively 
(Table 3) - are 2/arcmin^ and 0.7/arcmin^. Thus, given the opti- 
cal/160 yum ofl'sets (Fig. 1) the a priori probability of an acciden- 
tal alignment is very small, clearly smaller than 5%. We refer to 
a future paper (del Burgo et al., in prep.) for a detailed study on 
source contamination within the DUNES fields. 



4. Analysis 

The small PACS excesses above the photospheric fluxes are in- 
terpreted as due to cold debris discs around the stars. Black body 
dust temperatures, Tdust, can be estimated from the 100 and 160 
yum excess fluxes for HD 210277; in the case of a Men and HD 
88230, an upper limit for Tdust can be calculated taking the 100 
yum flux as 3cr statistical noise. Tdust for the three stars is < 22 
K (Table 4); the corresponding inner radii of discs with black 
body grains at this Tdust, considering the luminosity of the stars, 
are also given in Table 4. Fig. 1 shows the excellent agreement 
between the combined SED of the stellar photospheric fits plus 
22 K black bodies and the observed SEDs. 

Deconvolution of the images can be used to estimate the true 
size of the resolved discs in a Men and HD 88230. Our method 
first removes the photospheric contributions from each image 
by subtracting a PSF with a peak scaled to the predicted photo- 
spheric flux level. The PSF model uses a Bootis images, rotated 
to match the roll angle of the telescope during the observations of 
the DUNES stars. After star subtraction, the images are decon- 
volved using both modified Wiener and Richardson-Lucy algo- 
rithms. The noise model takes into account that the main contrib- 
utor is the telescope thermal emission. Both algorithms produce 
consistent results although with diff'erent noise patterns. Fig. 2 
shows the 160 yum star- subtracted and the Wiener-deconvolved 
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Table 3. Optical positions of the stars and of their PACS 100 identified counterparts. Observed PACS fluxes with Icr statistical 
errors (i^pAcs), and predicted photospheric fluxes (F*). Flux units are mJy. 



Star 


Optical position 


PACS 100 //m position 


Offset 




PACS lOOyL/m 




PACS 160yL/m 




(J2000.0) 


(J2000.0) 




^PACS 


/tlOO 


^PACS 




a Men 


06 10 14.47 -74 45 11.0 


06 10 14.53 -74 45 11.1 


0.3 


17.8+1.3 


18.0 -0.2 


14.4+2.0 


7.0 3.7 


HD 88230 


10 1122.14-^49 27 15.3 


10 11 21.88 +A9 27 17.2 


3.2 


22.5+0.9 


25.7 -3.6 


16.0+1.7 


10.0 3.5 


HD 210277 


22 09 29.87 -07 32 55.2 


22 09 30.08 -07 32 54.1 


3.3 


8.5+1.0 


4.6 3.9 


12.4+1.6 


1.8 6.6 




40 



-20 







40 



-20 





1 1 1 1 1 


HIP109378 


200 AU 

























40 



20 -20 

R. A. Offset (arcsec) 



-40 




Fig. 1. PACS images, isocontours and SEDs of a Men/HIP 29271 (left), HD 88230/HIP 49908 (middle) and HD 210277/HIP 109378 
(right). Images: 100 jim (left), 160 jim (right). North is up and East to the left. Isocontours: 100 jim contours are in green colour 
while the 160 yum ones are in red. a Men: 100 yum contours (10, 20, 40, 80 and 90 % of the flux peak); 160 jim contours (20,40,60,80, 
90 % of the flux peak) HD 882308: 100 jjm contours ( 10, 20,30,40,60, 80, 90 % of the flux peak); 160 jim contours (20, 30, 40, 
50, 60, 80, 90 % of the flux peak) HD 210277: 100 yum contours ( (30, 40, 50, 60, 70, 80, 90 % of the flux peak); 160 yum contours 
(50, 60, 70, 80, 90 % of the flux peak). The lowest contour in all cases is ^ 3cr. The optical position of the stars are indicated by the 
symbol in the isocontour plots; a segment indicates the projected linear sizes at the distance of each star. SED plots: black line 
is the photospheric fit while the blue line is the photosphere plus a 22 K black body. 



images of a Men and HD 88230; estimated sizes of the decon- 
volved sources are 16'' and 21'', respectively. Table 4 gives the 
linear sizes of the semi-major axes (~3cr contours) from both the 
original and deconvolved PACS 160 yum images. In the case of 
HD 210277, the value in Table 4 corresponds to an upper limit 
of 12", i.e., approximately the 160 yum beam size. The compar- 
ison of both the directly observed and deconvolved semi-major 
axes to the estimated radii from Tdust clearly indicates that the 
observed discs (assuming as a first approach symmetric discs) 
are smaller than the expected sizes for black body discs. 

Dust fractional luminosities, / (Table 4), can be obtained by 
taking the 160 yum fluxes and assuming the star temperatures 
given in Table 1 and Tdust = 22 K (Beichman et al., 2006). The 



estimated / values are of the order of 10"^, close to the Kuiper 
belt level. Essentially similar / values are obtained taking the 
maximum wavelength corresponding to 22 K, i.e., 230 yum, and 
its expected flux extrapolating from the one measured at 160 yum. 



5. Discussion 

The debris discs in this work are a new class of discs charac- 
terised by an excess at 160 yum, little to no excess at 100 yum, 
and no excess emission at shorter wavelengths. These discs are 
the coldest and least luminous ones known to date; they are sig- 
nificantly colder and fainter than other observed DUNES discs 
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Fig. 2. 160 jim star- subtracted (left) and deconvolved (right) im- 
ages. Up: a Men/HIP 29271. Contours: 10%, 20%,40%, 80% of 
the peak. Bottom: HD 88230/HIP 49908. Contours: 20%,40%, 
80% of the peak. North is up and East to the left. Eield size is 
60"x60". 

Table 4. Black body dust temperatures from the far-IR excesses, 
Tdust and the corresponding estimated radii. Observed Scr linear 
sizes (semi-major axes) of the extended emission at 160 fim in 
the original and deconvolved images. / is the dust fractional lu- 
minosity. 



Star 


^dust 


Size 


Size 


Size 


/ 






Est. 


Orig. 


Deconv. 






(K) 


(AU) 


(AU) 


(AU) 




a Men 


<22 


>147 


92 


81 


9.7x10-^ 


HD 88230 


<22 


>62 


56 


51 


1.6x10-6 


HD 210277 


22 


160 


<130 




5.4x10-6 



(e.g. Liseau et al. 2010, Marshall et al. 2011). Their interpreta- 
tion poses significant challenges. 

The shapes of the SEDs suggest that the dust is located in 
a ring with a larger inner void. We have seen that the observed 
disc radii are smaller than those implied from Tdust (Table 4). In 
other words, r^ust is smaller than the black body temperature, 
Tbb ~ 25-30 K, given the observed sizes and stellar luminosi- 
ties. Thus, we need to reconcile these seemingly contradictory 
results. Since one would reasonably expect a broad grain size 
distribution, which would include small grains hotter than those 
indicated by the estimated Tdust, the question is how to make the 
dust cold, while allowing a size distribution, but still keeping the 
disc radii within the observed ones. 

The low temperatures of grains require them to be large and 
highly reflective. The required albedo can be derived from Tdust 
and Tbb. Taking Tdust = 22 K as representative for the three ob- 
served discs, and Tbb = 25-30 K, the albedo would be >50%. Eor 
comparison, surfaces of trans-Neptunian objects contain signif- 
icant amounts of ice (Barucci et al. 2011) and many, especially 
large ones (e.g. Pluto, Eris), have albedos in excess of 50% (e.g. 
Vitense et al. 2010). Accordingly, Pluto has a surface brightness 
temperature ~ 10 K below the black body value (Gurwell et al. 
2010). It is natural to expect that dust released from the surfaces 
of such objects would have similar properties. 



In addition, one has to explain why small grains in the cold 
discs are depleted. One possibility is to assume a very low dy- 
namical excitation of dust-producing planetesimals, so that discs 
would be devoid of small particles (Thebault & Wu, 2008). The 
reason is that low collision velocities between large grains, un- 
afl'ected by radiation pressure, create an imbalance between the 
rates at which small grains are produced (low) and destroyed 
(high). Low collision velocities are compatible with low orbital 
velocities at the large radii of the cold discs. Besides, since the 
surface density of the solids far from the central stars is also 
low, planetesimal accretion scenarios predict very long accre- 
tion timescales (Kenyon & Bromley, 2008). Thus, large plan- 
etesimals that would excite the discs may have failed to grow. In 
fact, low dynamical excitation has been inferred for other large 
debris discs, e.g. HD 207129 (Lohne et al 2011). 

The tenet that the cold discs are probably in a low dynam- 
ical excitation state would be difficult to reconcile with the ex- 
istence of large planets in the discs, since they would stir the 
discs too strongly. The emission around a Men and HD 88230 is 
asymmetric, which might be suggesting the presence of a giant 
planet as in Fomalhaut (Kalas et al.2008). The problem can be 
mitigated if the planets are in nearly-circular orbits or the plan- 
etesimals have low eccentricities, as suggested for Fomalhaut 
(Chiang et al. 2009). HD 210277 hosts a planet at 1.1 AU with 
Cp = 0.47; however, based on the formulae by Mustill & Wyatt 
(2009), a stirring front from such a planet should not be able to 
reach the -130 AU-sized disc on Gyr timescales. 

Another point is the origin of the large inner voids in the cold 
discs. These could either be due to clearing by planets or may 
reflect the accretional and collisional history of primordial discs. 
In the latter case, the observed radii are those at which solids 
could reach "right" sizes, "right" degree of dynamical excitation, 
and/or were able to survive over the stellar age, to produce the 
observed emission. A detailed analysis of the nature of the cold 
discs is beyond the scope of this paper, and we defer to Krivov et 
al. (in prep.) where possible scenarios will be discussed in detail. 

6. Conclusions 

We have presented Herschel PACS observations of three stars of 
the OTKP DUNES sample. The observations reveal a new class 
of debris discs with fractional luminosities close to the Solar- 
System Kuiper's belt, but are colder and larger. These discs are 
a challenge to current models explaining debris discs around 
mature solar-type stars such as either the usual collisional- 
dominated disc scenario or low dynamical excitation discs. 
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